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Abstract 

This paper discusses a covariance study on the feasibility of using 

phase on short baselines to track the TDRSS satellites < 0rb ^ ccur ^ t s /^ n t ^ ^en 

station-differenced carrier phase data and range data collected from White Sands, NM are p c 
for various configuarations of ground stations and range data ^ 

nnsition accuracy of 25 meters can be acheived using two orthogonal baselines of 100 km 
So“Sl phase data and range data with 1 m accuracy. Relevant con^urauon 
parameters for the tracking system and important sources of error are examine . rdtS 

these data to redetemine the position after a station keeping maneuver is address^ inK BRTS 
sysmm Which is currently used for TDRSS orbit determination, is bnefly desenbed and its 

errors are given for comparison. 

I. Introduction 

SSumen^S 

reouired to be known to 50 m at one-sigma and their position must be redetermi e J 

involving GPS tracking techniques, are also being explored. 

In this paper, the possibility of using station-differenced earner ™ 
satellite to White Sands. The differenced phase between a pair of st t a “° c n f ^ S orccision of the 

differenced pha“ C obse^bk has a remaining bias due ^VpLn^ofrsky 6 poshio^oveTa data^amis 
through the two stations; as a result, only the change P interferometry (CEI) 3 in which the 
effectively measured. This technique differs from resolved. While station-differenced 

phase bias between the stations and the integer eye g y ODerat i ona ii y simpler and does not 

biased phase provides weaker position information than CEI, it.is operationally simp 

require a capability to record quasar signals for calibration as required for CEI. 
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— — — ruMuun nrror ( m) 

Computed error: " — 

Range data noise ^ 

10 m white noise, 10 m bias 


Doppler data noise ^ 

0.003 Hz white noise 


Total computed error ^g 

Considered errors: 

Solar pressure 

2% error in reflectivity 

4 

Troposphere 

5 cm zenith error 

2 

Ionosphere 

10 TEC error 

18 

BRTS station locations 

5 m error per axis 

13 

RSS position error: 

1" Thp rlfltrt haup O lO PAr* /"V rA : — 

73 


hour arc. 


Ini h i?i aper ’ b f? el i nes s ™ a11 enough to flt within the few hundred kilometer footprint of the TDRSS 
^ in ^ S j?i a ^ considered. Because of the high precision with which phase delay can be 

measured the station-differenced phase data type can provide g£od plane-of-sk7veK 
with baselines of this modest size. For stations within about 100 km, h 

drifts^n stpawte^ statilw 'docks^i 'Ba^elines°of Tkm W ' dl 

SSSS baSelineS 0f 100 km a " d 500 k ™ ta » hich the « sepS 

II. Current TDRSS Orbit Determination 

a) BRTS system 

7obe B FTrh C TnR S «S S ° f SeVCral ^“"d-based transponders at four near equatorial locations around the 

BRTS occupies the SA or MA user antenna of each satellite for about four minutes eve^ four HJ 

Because it is the operational data type used for TDRSS orbit determination, many studies have been done 
to determine the accuracy of the BRTS solutions.! BRTS regularly achieves sT mOm o^ siVrna 
posmon errors. For completeness, an error budget for BRTS is included here (Table 1) A data accuracy 
of 10 m for the BRTS range data and 0.003 Hz for the BRTS doppler dam is used for ^0 second 
integration times Four minute arcs of data are scheduled every four hours, and an epoch state solution 

a b?7 n7 3 3 t 4 h0UI k 313 arC ' As y stematic error in the ^ge measurements was included by estinL ing 
a bias parameter with an a priori error of 10 m. Error! from propagation media solar^ prasTrf 

ismo e ling and station location mismodelling are considered at the levels shown in Table 1 The 
computed error dominates the TDRSS position "error; however, ionosphe7mis m rel^ Hnl and s\aSn 
location uncertainty contribute significant consider error. g ana stall0n 

bj White Sands Ground Tracking System 

Sa 5? SltC | ha i s angle and range radio metric data available for rough TDRSS position 
on. e angle data, consisting of azimuth and elevation measurements from antenna pointing. 
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Table 2 

Differenced Troposphere Models for Several Baseline Lengths 


Station separation (km) Differenced troposphere Differenced troposphere 

r . , x i ah i a-1 r* An ctnnt I Cl 



have a precision of0 £ ^ ^ ^ ^ wLTfhe TORSS ^ ^skio^soludo'n's p^paga'ted 

?u?ent ranging system enhanced with station-differenced phase data is also studied. 

III. Data Modelling and Filter Assumptions 

at Configuration 

The Dosition solutions are found for 24 hour data arcs lnciuaing stai nA v T o f;i ter 7 Position 

with respect to White Sands. Hence, only the solutions for TDRS-east are show . 

ssress 

telt^ha^l7 , Wlh^of > a c^le^e ^ta we^fit t three^atio^^he C pliaM*bi'ase^ 

ambiguity and uncalibrated signal path delays at each station, are 
taken to be constant and initially unknown. 
h) Noise modelling for stati o n-differenced Phfli>$ dvUft 

ttvo lines of sigh, and the result mainly from 

fkS^o^^th^sil^l^^sation media. For 14 GHz signals, the dominant media fluctuations are 
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Clock 

Figure 1 . Configuration of antennae for carrier phase measurements used in covariance study. 



Clock type 


Rubidium 

Cesium 

Hydrogen Maser 


Table 3 

Station Clock Models Used 


Phase Random Walk Frequency Random 
growth (km/Vs) Walk growth 

((km/s)/Vs) 


Phase growth in 
24 hours 
(cm) 


1.5xl0’ 6 

1.5xl0' 6 

1.5X10' 8 


4.7xl0- 10 
1.5x10' 11 
4.7xl0* 12 


from troDOSDhere so charged particle media are ignored in this analysis. The data noise for differenced 
phase data from ’two statfonssharing a frequency reference is simulated by estimating a stochastic 
ffoposphere model at each station. To simulate the noise in the differenced phase data from stations with 
separate references, a stochastic clock model at each station is also estimated. 

il Troposphere modelling 

The sigma-squared and time constant of the Gauss-Markov stochastic troposphere model used for each 
baseline^are^ determined by fitting the autocorrelation of the station-differenced troposphere to an 
exponential (tabl^. In the cast of a 1 km separation between stations, two exponential . one 
accounting for the short term fluctuations and one for the longer term fluctuations, arerequired to fit t 

wen. For .he longer baseline,, a single exponential fits wel The — etaon 
function is calculated in a flat-earth frozen troposphere model using a wind velocity of 10 m/sec and a 
troDosDheric height of 1 km 8 ’ 9 The autocorrelation is calculated for north-south and east- west baselines 

for the 10° elevation of the TDRSS satellites from White Sands and the 
proiection factor between f he wind direction and baseline (fig 2). Since there is little dependence on 
baseline orientation, the autocorrelation with the largest value at t=0 ^ a given 

model all the baselines of that length. For symmetry, a troposphere model * aPP^« The Sffem^ed 
a sigma-squared half the value of the differenced troposphere autocorrelation sigma. The differenced 
troposphere fluctuations grow with baseline; but they grow slower than the enhanced precision in angle 
measurement due to the longer baseline. 

iil Station Clock Modelling 

Slightly longer baselines can be considered if frequency reference sharing between stations is not 
required. Baselines of 100 and 500 km in which the stations each have their own frequency reference are 
considered in this paper. If the stations forming a baseline have separate frequency references, drifts 
between the clocks at the stations increases the noise in the measurement of station-differenced phase. 

Typical frequency standard stabilities have a short term behavior of a white frequency 
term behavior of a random walk in frequency. In the OASIS program, the phase of a station clock is 
m^eUed as a polynomial in time, 0 = <Do + coo* +ax 2 - where x is time past some epoch^ 
parameter, coo is a drift parameter and a is a drift rate parameter In this analysis, die ^white 1 freqi ue y 
noise behavior of the station clocks is modelled by applying a random walk noise mode I ] t °^. t ^ e ^ 
bias and the random walk of frequency behavior is modelled by applying a random walk model to coo, 
the station clock drift Clock models representing the performance of rubidium, cesium and hydrogen 
mLer standards, showni^ table 3, are Ssed in L study.10 The rubidium and 

comparable short temt stability, but the cesium has better stability on a 24 hourorFeaKrtimesc^^A 
hydrogen maser is stable enough over a 24 hour data arc that would result tn position solutions 
comparable to the case in which the stations share a common clock. 
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Figure 2. Autocorrelation function of the troposphere differenced between two lines of sight with 

km^r^irn 1 ? 0 t 7 Station u s fo ™ in S an east-west or a north-south baseline of size: a) 1 km, b) 10 
km or c) 100 km. These are based on the model developed in ref 9. 
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Azimuth(°) 

Figure 3. Signature on the plane of the sky of the TDRS-east as seen from White Sands, NM. 


cl Range data 

Station-differenced phase data are able to trace the small signature in the plane-of-sky made by a 
geosynchronous satellite over a day (fig 3). In the geomedy at hand, this demnmnes 5 elen^ts of the 
TDRSS orbits well, but leaves the longitude-of-node poorly determined. In addition, an orbit solution 
with this data type alone is very sensitive to mismodelling of forces such as solar pressure. A sma 

mismodelling resets in a several hundred meter error in the satellite position s OTbiTdl ^temfnaS wUh 
satellite longitude. For these reasons, it is impractical to perform TDRSS orbit determination w t 

station-differenced phase data alone. 

In the analvsis ranee and doppler data from one station are included in the simulated data set. The range 
data^ ^are mSeied with a whhe noise measurement error along with systematic bias parameter which is 
estimated, with an a priori constraint corresponding to the ranging system .calibration accuracy^ Doppler 
data are modeled simply with white measurement errors. We will consider weights and biases for t 
range measurements between 1 and 30 meters. 

Range data can help determine the along-track position, which helps constrain the s ^ e ^ lon ^^ e e 
Low precision range data (20-100m) is good enough to control the error due to solar P res ^e 
mismodelling. Bette § r range data are required to reduce the computed along-track error to an acceptable 
level The range partial with respect to the longitude of the TDRSS satellite orbit, 8p/8<)>, for the 
geometry discussed here is 1/7, so a range measurement of better than 7m is required to get *e^ongi 
component of the TDRSS position error below 50m. The range data precision does not have tote quite 
this good, since it averages down over several measurments. However, the systematic error in the range 
measurment must be less than 1/7 of the required position error. 
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Table 4 

TDRS-east satellite position errors using station-differenced 
phase data from connected-clock stations andrange from White Sands. 


inputs 

Data weights: 

Station-differenced phase baseline (km) 1 * * * " 
Dopper Noise (mm/s)* 

Range Noise (m)^ 

Range Bias (m) 

Consider Inputs: 

Range station location(m) 

Station-differenced phase station locations(m) 

Solar pressure(% of reflectivity) 

results (epoch) ” 

Computed error (m) 

Consider errors (m) 

Solar pressure 
Range station location 
Station-differenced phase station location 
RSS position error (m) 

t Two orthogonal baselines of this size are used, 
t Data noise for 10 minute points. 


case 1 

case 2 

case 3 

case 4 

case 5 

100 

100 

100 

10 

i 

1.0 

1.0 

1.0 

1.0 

1.0 

10 

10 

1 

1 

1 

30 

10 

1 

1 

1 

2 

2 

2 

2 

2 

.35 

.35 

.35 

.35 

.35 

2% 

2% 

2% 

2% 

2% 

247 

84 

14 

25 

63 

12.1 

10.4 

12.4 

12.4 

12.6 

14 

14 

14 

14 

14 

0.12 

0.12 

0.10 

0.12 

11.9 

248 

86 

23 

31 

67 


IV. Covariance Results 

a) Computed Errors 


Inhn S 1 *! e COmpilt f d e 'T or ? ! hat rcsult in several scenarios with station-differenced phase on two 

i) Dependence on Range Bias 

As is apparent from Table 4, the position accuracy is limited by the ranging system accuracy in several 

the 'n g H r f 4 ShOWS t! he dependence ofTDR S position error on range bits a priori un^enainty when 

the range data are combined with stanon-differenced phase data from various sized baselines. In ail cases 

' h i ™, ge t data noise ^ as 1 JP - In order tha t the position error be limited only by the station-differenced 
phase data, the range data bias must be better than approximately 1 m in the 100 km baseline connected 
frequency reference case. For smaller baselines, the station-differenced phase data are weaker and the 
TDRS position error becomes dominated by the station-differenced phase data at larger values of range 
0 1 casc £, ran g e data largely determine the component of the TDRS position in the 
to rfh f tUde direcDo /V whll< ; ^ station-differenced phase data determine the other orbital elements. In order 
to THR st T n§e b n S ° f m U may ^ necessar y to upgrade the current 2- way ranging from White Sands 
«v«rem fh ypiCa y ’ a ran 2 in S system can be calibrated to about 90% of the inverse bandwidth of the 
system, thus a ranging system calibrated to 1 m would require a 30 MHz bandwidth. 
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Figure 4. R.S.S. position error for a TDRSS satellite using a one day arc of station-differenced ptra and r ^8 e da ^ fr0I J 
White Sands as a function of systematic error in the range measurement. The range noise is held fixed at 1 m 

minute point 


til Dependence on Troposphere model 

Since the troposphere model is the important source of data noise for the station-differenced phase damn 
the C case where the station s share a common clock, it is useful to know how the troposphere model 
Darameters affect the TDRS position error. The two parameters that describe the station-differenced 
troposphere as modelled in this paper, are the steady-state sigma and the correlation time constant. 
Figured shows how TDRS position accuracy varies with these model parameters. In Fig. 5a, The 
XU= o°f*e troposphere sigma is varied by a factor of two from the nomtntd v^es shown tn Table 
2 while the time constant is held fixed at its nominal value. In Fig. 5b, the time constant tor tne 
stochastic station-differenced troposphere is varied, while the steady-state sigma is held fixed. In bot 
cases the position error changes slowly with the variation in troposphere parameters. 


jin Effect of S eparate Station Clocks 

Because a configuration in which the stations have independent freqeuncy references may be j® 
S we consider it as well, despite the cost in position accuracy. Position determination with stauon- 
differe^nced^ohas^ depencte on tracing the signature in the plane-of-sky from the station made by the 
“like In the c“se Ta geosynchronous sarfllite like TDRS. there is a 24-hour period to the s tgnature 
(Fig 3), so that is the time scale during which the difference in the station clocks .must be t stable The 
expected plane-of-sky error from clock phase error growth is approximately pA<)>c/B, where p 
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Figure 5 a.) R.S.S. position error for a TDRSS satellite using a one-day arc of differenced carrier-phase 
and range data from White Sands as a function of sigma of station-differenced zenith troposphere, b.) 
R.S.S posmon error for a TDRSS satellite using a one day arc of differenced carrier-phase and range data 
from White Sands as a function of station-differenced zenith troposphere autocorrelation time constant. A 
range data noise of 1 m and range bias of 1 m were used in both figures. 
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Figure 6 shows the expected position error as a funC "^ 

measuring phase (see Fig 1) have then own each station has a 

^""hlSteSge data have a 1 m data noise and 1 m bias error and 
are unaffected by the clock model. 


fO Consider Errors 

A consider analysts was performed to assess the aensiuvuyoj 

and no, included in Table 4. The same consider errors are 
applied to all the cases described in this paper. 

In this covariance analyis ^lar^ressure is sofa^fl a x!°rTs the d Sa ween^h e sun and the 

satellite surface. F — CA (1 + tVk* refVctine area is assumed to be constant 

wifh^vah/e of^ /o^^^Sol^pressure^^smodelling is simulated as a 2% of the solar reflectivity 
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lputs 

carrier phase baselines (km) 
range weight (m) 
range bias (m) 
doppler weight (mm/s) 
results 

error after 30 mins^ 
error after 60 mins^ 
error after 90 minst 
error after 120 mins* 


100 

1 

1 

1 


10 

1 

1 

1 


379 

249 

248 

247 


301 

94 

91 

88 


287 

124 

319 

49 

97 

242 

43 

85 

206 

40 

76 

182 


o ^ ui Udid irum oeiore me maneuver and data after the man 

column one of table (see text). The errors include consider errors at the same levels as table 2. 

Wlti L a Homina 1 v al ueof 1.42. It has been shown that solar pressure effects on the 
™ SS s ^ ellltes can be modelled to 2%M Operationally, the solar reflectivity might be estimated 
enlarging the computed error somewhat, and a much smaller mismodelling error could be considered. 

The error due to uncertainty in the antenna locations is also determined with consider analysis A 2 m 
uncertainty in each component of the ranging station location is used. A 35 cm uncertainty in each 
“’ npone f of l J e stations measuring carrier phase is used. This corresponds to allowing about 0 5m 
uncertainty in the baselines formed by pairs of stations. The sensitivity^ to the baseline length error is 
small except in the case of the smallest baselines studied. e gtn error 1S 

c) Recov ering Position After a Manpnv^r 

The TDRSS satellites occasionally make corrective maneuvers in order to stay in their desired orbits to 

tlT era " C n The rei ’ uir f mem for TDRSS is 50 "> onLigma whhi^wo houI? 

^ , h ,"l aneuver - Because 8°«J position determination with station-differenced phase data reauires 
£>"' “ S “ ‘ s ™P ossib ^ to redetermine all the state parameters to tw“aSc? 3 ,wo 
hours. Instead of redeteimtntng the whole state, we may take advantage of the good instantaneous Diane- 
of-sky velocity information in this data type by estimating the velocity change in the satellite state 
associated with the maneuver. Line-of-sight doppler data collected along with the ™Ye dam at one 
station, provides the third component of the satellite velocity. 

It is assumed that there is nearly a day long arc of station-differenced phase, range and doppler data prior 
e maneuver. The corrective maneuver is modelled as a velocity impulse with well known time of 
bum. No a priori knowledge of the error in the impulses is assumed. Table 5 shows posTrion^mor 
. ieva . be after such a ™ ane uver, if the state and the maneuver are both estimated using a day long arc of 
data ending wo hours after the maneuver. In each case, the position error is about twife as l^rge Ste 
a full day data arc unmterupted by a maneuver. s 

V. Conclusions 

We have examined the possibility of using station-differenced earner phase with stations forming small 
?0 and7^ satellites. We find that fwo orthogonal baselines Sen 

^"-differenced phase data alone poorly determine the 
ongitude of the spacecraft, and thus range data with systematic error of about 1 m must be included to 
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baselines are required. 

The station-differenced phase data along with two-way doppler ran be used “ 

change associated with a s ! atl ^ n ' ke ^ to recover the best possible position 

accurac^'with this meS No"" towfedge of Je size’of the impulse is required; though the time 
of the impulse is assumed to be known. 

Range data with accuracy of about 1 m is crucial 

the exact vaLe of the 

bias. 
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